Effects of hadronic loops on the direct CP violation of B c 
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It is well known that the final state interaction plays an important role in the decays of B-meson. 
■ The contribution of the final state interaction which is supposed to be long-distance effects, to 

the concerned processes can interfere with that of the short- distance effects produced via the tree 
and/or loop diagrams at quark-gluon level. The interference may provide a source for the direct CP 
violation Acp in the process — > D°tt + . We find that a typical value of Acp when the final state 
interaction effect is taken into account can be about —22% which is different from that without the 
final state interaction effect. Therefore, when we extract information on CP violation from the data 
which will be available at LHCb and the new experiments in B-factories, the contribution from the 
final state interaction must be included. This study may be crucial for searching new physics in the 
future. 
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I. INTRODUCTION 



One of the most intriguing goals in the high energy physics is to look for new physics beyond the Standard Model 
(SM) via heavy hadron production and decay processes. The reason is that new physics which generally has a higher 
energy scale may be observed at the processes involving heavy flavors. Among all the possible quantities which are 
^\ ' experimentally measurable, CP violation provides a more sensitive window to the new physics effects. Direct CP 
t— I , violation at i?-physics has been observed by the Babar and Belle collaborations [3, Q , which is indeed a great success 
after confirmation of non-zero e'/e at K-systems. Another promising place to study CP violation is the meson _B C , 
which is composed of different heavy flavors. 
£ — ■ Since the CDF Collaboration observed B c meson in the semileptonic decay B c — > J/tp + I + v studies on B c 
' have drawn great interests from both theorists and experimentalists of high energy physics. Decays of B c can be 
realized via b— decay, c— decay and annihilation of b and c 0|- Many theoretical works have been dedicated to study 
the decays of B c [1, H, 0, [1| ■ A relatively complete discussion about its spectrum, production and decays was presented 
in a review Q. Because of the specific characteristics of its decay modes, the direct CP violation is an important 
observable which may provide valuable information towards the mechanism governing the transition and probably 
unveils a trace to the new physics beyond the SM. 

In this work, we are just looking for a new source for the direct CP violation in Bc decays. The direct CP violation 
is caused in general, by an interference among at least two channels which have the same final state, but different 
weak and strong phases. The CP quantity Aqp is proportional to 

. 2|AiP 2 |sin(0i -0 2 )sin(ai - a 2 ) 

Acp ~ 



|Ai| 2 + \A 2 \ 2 + 2\A 1 \A\ 2 \cob(0i - 2 ) cos(ai - a 2 ) ' 

where A\,A 2 are the amplitudes of the two distinct channels and 0i,0 2 , a\,a 2 are their strong phases and weak 
phases respectively. 

These phase differences coming from either quark level or hadron level. At the quark level the strong phase difference 
usually occurs via the absorptive part of the loops involved in the calculation. The strong phase may also occur at the 
hadron level. As a matter of fact, it is well known in the kaon system. When one studies the direct CP violation, i.e. 
e'/e, the phase shifts in the 7T7t scattering provide the strong phase which are necessary to result in CP violation. But 
recently most of the works to study direct CP violation concentrate on the strong phase induced by the absorptive 
part of the loops. Especially, the strong phase is coming from the absorptive part of the penguin diagram(s) which 
contribute along with the tree diagram to the amplitude. In that case the CP violation is induced by the interference 
between the contribution of the tree diagram and that of penguin. 
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The total width is related to \A\ + A-^ . If one of the amplitudes is much smaller than the other one, the width 
should be only depend on the larger one, say A\, thus one can ignore the smaller one when he is calculating the decay 
width. However, even though \A-x\ <C \A\\, the numerator of Aqp is proportional to their product, so one cannot 
ignore the smaller contribution, otherwise he would get null CP asymmetry. 

In that case, obviously the contribution from the penguin diagram is much smaller than that from the tree diagram, 
so that if only the decay width is needed, one can completely ignore the contribution of the penguin. However, for 
evaluating the CP violation, he by no means can dismiss the penguin contribution. 

In the SM, the weak phase originates from the Cabibbo-Kabayashi-Maskawa matrix, and the strong phase is induced 
by the absorptive part of loops. At the quark-gluon level which is responsible for the short distance effects, the strong 
phases may originate from the absorptive part of loops, for example, the penguin diagrams. On the other aspect, 
the final state interaction (FSI) plays an important role in B-physics, as fully discussed in literature [lol ]. At the 
short distance, the direct CP violation usually is caused by an interference between the tree-level contribution and 
the loop-induced one because they have different weak and strong phases (in fact the tree diagrams do not contribute 
a strong phase). Therefore an interference of the long-distance contribution with the short-distance ones may change 
the theoretical prediction on the CP violation. In fact, the FSI effect is extensively applied to the discussion of the 
CP violation of B and D decays [lU, E3, ■ 

Indeed, by the quantum field theory, the lagrangian can be a combination of various pieces and each of them 
corresponds to different processes. For our transition matrix element M —out (f\i)in, one has 

M = (f\LpQ CD + T[LhadL^p QCD ]\i), 

where LpQQ D corresponds to the lagrangian which includes QCD and weak or electromagnetic interactions, the 
superscripts (1) and (2) denote the lagrangians which can lead to different final states, whereas L^ a d is the lagrangian 
at hadron level. Then we further write the matrix element as 

M = (f\ L PQCD\i) + ^{f\Lhad\n){n\L% CD \i), (1) 

n 

where the intermediate states \n) are a complete set of hadrons with proper quantum numbers and the matrix element 
{f\Lhad\n) is just the hadronic scattering process and corresponds to the hadronic loops in our work. 

Generally, the long-distance effects due to the FSI refer to the rc-scattcring of the intermediate hadrons which 

emerge at the direct decays, into the concerned final state and it is depicted by the term (n\L P Q CD \i) . In these 
channels with the intermediate hadronic intermediate states may have different weak phases from that of the short- 
distance production channel occurring at quark-gluon-level. In the re-scattering processes, phase shifts exist due to 
strong interaction and thus can offer strong phases. And an extra strong phase which is definitely different from that 
induced by the quark-level loops, occurs from the hadron re-scattering processes (f\Lh a d\n). 

Since the loop contribution is suppressed by the loop integration, generally the second term of the above equation 
is smaller than the first one which we may call it as the "tree" level contribution( but maybe not the tree diagram in 
the common sense). 

The traditional PQCD calculation only takes care of the first term (f\LpQ CD \i) and (n\LpQ CD \i) in the second one, 
but leaves the part (f\Lh a d\n) to be dealt with in other theories, for example, the chiral lagrangian and etc. at the 
hadron level. This picture is clearly depicted in Cheng's paper It indicates that unless the "tree" contribution 
(i.e. the first term) is suppressed by some mechanism, the first term corresponds to the direct process, so that is 
always dominating for the total amplitude. If we only need to consider the total decay width, the second term may 
contribute a smaller portion (sometimes it might be enhanced by some mechanism, but generally is much smaller). 
However, as we deal with the CP violation and need at least two different channels, their interference enforces us not 
to abandon this term even though it might be much smaller than the first one. Indeed, one may argue that the loop 
diagram, such as penguin can also contribute a strong phase and interfere with the tree contribution to result in a 
direct CP violation, the loop contribution may have a similar order as we considered here and possibly even smaller. 
At least as we state above, we are looking for a possible source of CP violation in Be decays, i.e. the hadronic loops 
may contribute strong phases and cause sizable effects on CP violation as our numerical results given in the paper 
indicate. 

Therefore we would say that the PQCD framework works well, but we instead are looking for a supposed-to-be 
smaller effect which can result in observable CP violation. If there is a small double-counting (could be), that is 
because the wavefunction adopted in the calculation is not well defined. In fact because \n) generally are not the 
same as |/), the double-counting does not appear. 

At present the direct CP violation in B c decays due to short distance contribution has been studied by many 
authors [TJ, EH, EH, E3, E3 • But so far, the studies of the FSI effects on the direct CP Violation of B c are absent. In 
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this work, by taking into account the long distance contribution caused by the FSI, we would re-evaluate the direct 
CP violation in the decays of B c . Namely, we add a new contribution to the amplitude which has different strong and 
weak phases from that of short distance contributions which were calculated by many authors, thus their interference 
will significantly change the value of CP asymmetry in decays of B c . 

Indeed, before one can claim a discovery of new physics, he must exhaust all possibilities which the SM can provide. 
Therefore this work is also serving for the purpose and to determine if the FSI can result in a sizable contribution to 
the direct CP violation of B c . 

We choose the channel B+ — > D°tt + which should be one of the dominant decay modes of B c . In this channel, 
there exist hadronic intermediate states which are mainly composed of and J/tp. 

This paper is organized as follow. We present the formulation about 73+ — > D^*' + J/^ — > D°n + in [Til Then we 
present our numerical results. The last section is a short conclusion and discussion. 



II. FORMULATION 



The effective Hamiltonian related to B c decays is [1 

G 



ff = ^^^{c i (/ J )(c6)y_A)(dc)y-A+C 2 b ( M )(J6)y_A)(cc)v-A} 

+ ^V ub vl d [c b M{ub) V -A){du) v „ A + C b 2 {ri{db) v _ A ){uu)v- A }, 



(2) 



where the subscript V ~ A denotes the left-chiral current 7 M (1 — 7 5 ). C\ 2 (/u) denote the Wilson coefficients . 

Firstly we calculate the transition amplitude of B c — > J /ip at the quark level and the hadronization would be 

described by a few phenomenological parameters. The definitions of the relevant hadronic matrix elements are 

(0|^|P(fc)) = -ifrkp, (3) 

(OI^IVCM) = fv^mv, (4) 

where f-p and /y respectively stand for leptonic decay constants of pscudoscalar and vector mesons, fc^ is the four- 
momentum of the concerned hadron and e M denotes the polarization of the vector meson. One has = 317^(1— 75)52- 
In addition, the hadronic matrix elements of B c transiting into two mesons can be expressed in terms of a few form 
factors as 

(V(k 2 )\J^\B c (k 1 )) = P»MQ 2 ) + Q»f-(Q 2 ), 



-^{V{k 2 ,e)\J ll \B c {k 1 )) 



J {is^P a Q p F v {Q 2 ) - gT{P ■ Q)F A (Q 2 

mi + m 2 

+ pnp»F A {Q 2 ) + Q"P U F A (Q 2 )} 



with P M = (fci + k%) t : l and — (fci — kz)^- With the above formulas, we obtain 



(5) 
(6) 



and 



M[B+(p) - D+(pi) J/V(p 2 )] 
V2 ^ m Bc + m 



9°\V + P2) ■ (p - P2)F A (q 2 1 ) + (p + P2r(p + P2) A F^( q i) 



+ (p - P2Y{p + P2) X F A {q 2 )\ + a 2 f m [(p + Pi) X f + (q 2 2 ) + (p - Pi) X f- (<? 2 2 )] }, 



M[B+(p) -> D*+(p!) J/V(pa)] 



VcbV r 



D *m D , 



is 



rriB c 

?Ai „ 2 



\p+P2)r{p-p2)sF V (q 2 ) 



iG F 

ST(P + P2) ■ (P - P2)Ff(q() + (p + P2 r(p + P2YF*(qt) 



+ (p-P2Y{p + P2YF A (q 2 ) +a 2 /^ 



IS 



\p + Pi) rip - Pi)sF v (ql) 



m Bc + m D 

-g^iP + Pi) ■ (p - Pi)F A (q 2 2) + (p + piY(p + piYF+(ql) + (p - piY(p + piYF A (qI) 



(7) 



(8) 



with q\ = p — p 2 and q 2 = p — p\. The values of ai i2 will be given in next subsection. 
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A. Absorptive part of hadronic loop for J3+ — » _£)(*>+ J/tp — > _D°7r+ 

Now let us turn to evaluate the contribution from the long-distance effects which occur at the hadron level. The 
diagrams shown in Fig. [T]depict sequent processes — > £)(*)+ J/ij) — > D°7r + . 




FIG. 1: The final state interaction contributions to — » D°7r + . 
The effective Lagrangian at the hadronic level is suggested to be in the following forms as [2C 

Cd'w* = -gvi>.«ei a " x f > dj>tird a fy, 
C^vv = ig^-D^^VV -V&"D), 



(9) 
(10) 
(11) 

(12) 



with 7T = t ■ 7T, where fields X>M and X>W arc defined as X>M = (£)W°, DW+) and £>W T = (ijW , flW"). 

The process shown in Fig. 1 (a) is — > D + {j>i)J/ij)(p 2 ) — > 7r + (p3)D°(p4) where D*° is exchanged at t-channcl, 
and its amplitude reads 



Abs^ = - 



{2Kf5\m Bc - Pl -p 2 )M[B+(p) -> D + (pi)J/V(P2)] 



2 y (2^)32^1 (2tt) 3 2E 2 

x [ - flBT^ipl] [ - ig.j/^D'D^ al3 {.-ip2n)iq a ] {-g\v + P2xP 2 2u ) 



(13) 



Obviously, the conservation of angular momentum demands the contribution from Fig. 1 (a) to be zero. 

The amplitude corresponding to the process of — * D* + (pi)J/ip(p 2 ) — ► n + (p-i)D^(pi) where D° is exchanged at 
t-channel reads as 



Abs^ 



{2ir) A 5 /l (m Bc ~p 1 -p 2 )MlB+{p) J/^fe)] 



1 r d 3 px d 3 p 2 

2 J {2vf2E 1 {2ir) 3 2E 2 

x [ - gD>Dn(-ipt)] [ - g^DD{ivi - kT] (-g<r$ + V1<T £ H ) 

x(-3o)/iH — )-n 11 i m D ■ 

m tf> 1 ~ m D 



(14) 



For Fig. 1 (c), — > D* + (pi)J/ip(p 2 ) — ► 7r + (p3)Z?°(p4) where Z?*° is exchanged at t-channel, the amplitude is 
1 A c? 3 pi d 3 p 2 



-(27t) 4 J 4 (to Bc - Pi -U2)M[B+(p) D* + ( Pl )J/iP(p 2 )] 



2 J (2^)32^1 (2tt) 3 2£;2 
x [-ig D *D^e^ at3 (~ipi^)(-iq a )} [ - i3.//^D*_D£ ?AK ''(-«P2e)^«] 



, . <?/3g P w . Pi^ii/w . P2up 2 \, 
X\-90p H 2^)(-fl<7v 4 2 — H-ffwA -I — J 



(15) 
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In the above amplitudes, q — P3 — pi, and J-(q 2 ,rrii) etc. denote the form factors which compensate the off-shell 
effects of mesons at the effective vertices and may be described by the possible pole structures flfj | 

/ A 2 2 \ n 

^ q2 ^={l^t) > ( 16 ) 

where A is a phenomenological parameter to be determined. As q 2 — > the form factor becomes a number. If A 3> m%, 
it becomes a unity. As q 2 — > oo, the form factor approaches to zero. It reflects the fact that as the distance between the 
mesons becomes very small, their inner structures would overlap and the whole picture of hadron interaction breaks 
down. Hence the form factor vanishes at large q 2 and effectively plays a role to cut off the ultraviolet divergence. The 
expression of A is suggested to be [l(j 

A(mi) = m, + aA QCD , (17) 

where m,, denotes the mass of the exchanged meson and a is a phenomenological parameter. In this work, we adopt 
the dipole form factor T(q 2 ,mi) = (A 2 - m 2 ) 2 /(A 2 — q 2 ) 2 . 



B. The dispersive part of B+ -» D { * )+ J/ip -> D°n + 

In the above subsection, the absorptive part of the triangle diagram to the amplitude of the sequent process 
7r + can be easily obtained from the integrals (fT5|) , (|14[) and (fT5|) . The dispersive part of 
7r + can be related to the absorptive part via the dispersive relation [n], Hl[ 

Dis[B+ - Z?°tt+] = - / 7 a [ ds. (18) 

However, the cutoff which is phenomcnologically introduced and the complicated integral in eq. (|18p would cause 
unavoidable uncertainties to the dispersive part. In some of the former works, for estimating the decay width, the 
contribution of dispersive part was assumed to be small comparing with that of the absorptive part and ignored. 
However, for the direct CP violation, we must estimate the dispersive part and determine the strong phase induced 
by the triangle diagram, otherwise the strong phase would be exactly n/2. 

In this work, adopting the method in our previous work [22| . we obtain the dispersive part of — > !)(*)+ J/tj) —> 
D°ir + by directly calculating the triangle where the intermediate hadrons are not on their mass shells. The amplitudes 
corresponding to the process of £>+ — > D* + (pi)J/ifj(p2) — ► tt + (ps) D° (p±) where D° or D*° are exchanged are 

Dis^ = J ^M[Bt(p) ^ D* + ( Pl )J/^( P2 )][ - g D , D A~i P l)}[ - g^DDiiPi - iqr] 



x(-fttf)(-3<-/*)-3 2 T^ 2 b 2 ' m c], (19) 

Pi - m D P2 - m J/4 , q z - m z D 



and 



Dts(C) = J ^MBHP) ^ D*+(p 1 )J/i;(p 2 )}[-tg D , D ^e^(- i p^)(- l q a )} 

[ ~ ig.j/^D'D£ iXKp (-'ip2^iq K ] (-fl , /3p)(-fl , ^)(-5uA)- 2 — " 

p(-m D p 2 - m J/4 , 

x - 2 J*[cf,m 2 D .]. (20) 
q z — mi-). 

Due to the existence of the dipole form factors T 2 \q 2 , m 2 D \ and T 2 \q 2 , m 2 D ,] the ultraviolet behavior of the triangle 
loop integration is benign. These form factors play an equivalent role to the A— related terms introduced in the 
Pauli- Villas renormalization scheme [23|, [24[. Because the final expressions of eqs. (fl9|) and (f20"|) are complicated, we 
would collect some useful formulas in appendix. 



C. Direct CP violation 



The observable direct CP violation is defined as 

A C p = 



\M\ 2 -\M\ 2 
\M\ 2 + \M\ 2 



(21) 
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with 

M = M dir {B+ -> D\+) + M FSI (B + -> £> tt + ), 
AT = M dir (B- -> L)°7r-) + X FS/ (S C - -> D°7r-), 

where AA dtr (B+ — > D°7r + ) and A4 (B~ — > L>°7r~) were calculated in the approach of PQCD by the many authors 
[lH and they are written as 



M d "'(B+ -> D°n+) = V U {T U + P)[l - ze l{ -~< +S) ] 
M dir (B- -» 5°7T-) = F M *(T„ + P)[l - 2e i (T+ 5 )] 



(22) 
(23) 



with 





T c + P 




T u + P 



and 5 = arg 



T r . + P 



T u + P 



where V u — V u dY* b are the Cabibbo-Kabayashi-Maskawa entries, |T4/V^| 
7, A, z and S are given in Ref. [l8j . which are listed in Table U 



13^72 I Kb/Kb | ■ The values of T u , c , P, 
The amplitude of i?+ — + _D°7r + induced by the FSI effect which is denoted by the subscript " FSI" is: 



M FSI (B+ -> D\+) = Dis + i J2 Abt 



(24) 



j—a,b,c 



and 



M FSI (B+ DV 



M FSI (B~ 



D°%-). 



(25) 



T u 


22.621 +0.863i 


S 


123° 


T c 


-0.83 + 3.57i 


7 


55° 


P 


-0.474 - 1.722i 




0.085 


z 


0.28 







TABLE I: These values are taken from Ref. [T^]- Here T UtC and P are in unit of 10 -3 GeV. In this work, we need multiply a 
factor \J rrv'g^GF / \/2|k| ~ 4.83 x 10~ 4 to T u , c and P, because the formula for the decay widths adopted in this work takes a 
different normalization from that in Ref. [l8l ]. 



III. NUMERICAL RESULTS 



The input parameter set which we are going to use in this work, includes: bib c = 6.286 GeV, mji^ = 3.097 GeV, 
m D+ = 1.869 GeV, m D ,+ = 2.01 GeV, m D o = 1.865 GeV [H; U = 405 ± 17 MeV [H; Id = 222.6 d= 16.7±1;| MeV > 
f D . = 245 ± 20± 3 2 MeV j^. Vud = 0.974, V cd = 0.230, V cb = 0.0416 g D . D „ = 17.3, g D , D ^ = 8.9 GeV" 1 , 
.9 d Dtp = 7.9, g D , DlP = 4.2 GeV" 1 j27|; oi = 1.14, a 2 = -0.20 0. V ub = AX 3 (p~n]) = 0.00218 - 0.00335i. The 
wolfenstein parameters of CKM matrix elements: A = 0.2272, A = 0.818, p = 0.221 and fj = 0.340 with p = p(l - §) 
and fj = 7/(1 - §). Gf = 1-16637 x 10~ 5 GeV" 2 

The form factors in processes -B c — * and P c — * J/ip possess pole structures [1,0] 

F ^ = i "77772 ( 26 ) 

1 — a£ + oi, 

with C = q 2 /m 2 B . The values of -F(O), a and 6 are evaluated by some authors and for readers' convenience, we list 
their results in Tabic HT1 
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/+ /- 


F$ Fd 4 Fv 


D 


F(0) 


0.189 -0.194 


- 


a 


2.47 2.43 


- 


b 


1.62 1.54 


- 


D* 


F(0) 


- 


0.158 -0.328 0.284 0.296 


a 


- 


2.15 2.40 1.30 2.40 


6 




1.15 1.51 0.15 1.49 


J/i> 


F(0) 




0.66 -1.13 0.68 0.96 


a 




1.13 1.23 0.59 1.24 


6 




-0.067 0.006 -0.483 -0.002 



TABLE II: The values of F(0), a and b in the form factors of B c -> L> ( *> and B c -> J/tp [10]. 



20 




O 

-20 - 



-30 - - 

.40 I 1 1 . 1 . 1 . 1 1 

0.5 1.0 1.5 2.0 2.5 3.0 

a 



FIG. 2: The dispersive part and absorptive part of the amplitudes of — > D^*' + J/ip —* D°tt 



In Fig. [21 we plot the dispersive part and absorptive part of the amplitudes of — > D^ + J/ip — > D°tt + versus 
a which is allowed to vary within a = 0.5 ~ 3. In Fig. [31 we also list Acp with several typical values of a. For a 
clear comparison, in this figure, we also give the value of Acp calculated in Ref. [llj] by the PQCD approach which is 
purely induced by the short distance contribution (without considering the FSI). For clarity we also list some typical 
values of Acp with various a in Table III. 

IV. DISCUSSION AND CONCLUSION 

Recently the direct CP violation in .B-decays has been observed and it is expected to open a window for exploring 
new physics beyond the SM by which all theorists and experimentalists feel very inspired. Obviously, investigation of 
direct CP violation at B c decays would be of special interests because it is composed of two heavy flavors and may be 



a 


0.5 


1.0 


1.5 


2.0 


2.5 


3.0 


Acp 


-30.2% 


-28.1% 


-25.3% 


-22.4% 


-19.8% 


-17.5% 



TABLE III: The typical values of Acp 
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-16 




-32 I ' 1 ' 1 1 1 ' 1 ' 1 

0.5 1.0 1.5 2.0 2.5 3.0 

a 



FIG. 3: The direct CP violation. Solid line and Dash line correspond to the direct CP with FSI effect and without FSI effect 
respectively. 

more sensitive to new physics. On other aspect, before one can claim to find a trace of new physics, he must exhaust 
all possibilities in the framework of the SM. As indicated in literature, the FSI play an important role in B decays, 
therefore one has a full reason to expect that it is also significant at Be decays. In this work, we carefully study the 
contribution of the FSI to the direct CP violation via its interference with the contribution from the short-distance 
effects which arc induced by the tree and loop diagrams. Concretely, in this work, we calculate the amplitudes for 
_B+ — > D°tt + via sequent processes — > £)(*)+ J/tp — > D°n + and determine its strong and weak phases. 

Here we need to add some interpretation about application of PQCD. Even though we indicate the significance of 
the FSI for evaluating direct CP violation in Be decays, their absolute contribution is much smaller than that from 
the direct process which is calculated in the framework of PQCD. Therefore if only the decay width of Be is needed, 
one can ignore the contribution from the hadronic re-scattering, but as the CP violation is concerned as we see above, 
its contribution might be significant. 

Our numerical results indicate that the typical value of Acp with FSI effect is about —22%, which is different 
from the value —30.7% estimated in the PQCD approach without FSI [HI. On other aspect, one can also observe 
from Fig. 3 that the effect of the hadronic re-scattering on Aqp may change quite diversely depending on the input 
parameter. Especially, as one adopts a = 0.5, Acp is about —30.2% which only slightly deviates from the value 
obtained in the framework of PQCD, however, as a = 3, (even though a — 3 seems too large to be very reasonable, 
this effective coupling indeed can exceed 1 for hadron interaction, in fact, in some applications its value is set to be 
very large for fitting data) Acp would change to —17.5% obviously deviates from the value of PQCD. It indicates 
that the contribution of FSI to Acp is of opposite sign with that from the quark loops and the cancellation may cause 
remarkable effects when one analyzes the data achieved in a rather precise measurement. Thus our conclusion is that 
the contribution from the FSI is not negligible. 

In the future experiments, especially the LHCb, a great amount of data on B c will be accumulated and one may have 
a possibility to measure the direct CP violation of B c . If non-zero Acp (almost definitely yes) is well measured, one 
can look for a trace of new physics by comparing the measured value with the theoretical result. When one compares 
the data which will be available at LHCb and/or other experiments with theoretical predictions, the contribution 
from the FSI must be included. This observation may be crucial for searching new physics in the future. 
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Appendix 



Some useful formulas in the calculation of eqs. (fT9j) and ([20 



Here 



and 



i 



(2tt) 4 {p\ - m?)(p2 - m 2 )(q 2 - to 2 ) V g 2 - A 2 



A 2 - m 2 \ 4 



(4tt) 2 



(A 2 - TO 2 )?/ 1 1 

A 2 (toi,TO2,A) A(toi,to,2,A) A(mi,m2,TO) 

6\„,3 , 



(-A 4 - m 4 + 2m 2 A 2 )y 2 (A 6 - 3to 2 A 4 + 3to 4 A 2 - m 6 )y 



A 3 (toi, m 2 , A) 



A 4 (toi,to 2 , A) 



d 4 q I 2 /A 2 — TO' 

(27r) 4 (p 2 - m\)(p 2 - m 2 )(q 2 - to 2 ) V g 2 - A 2 

2j/(A 2 -to 2 ) rA(TOi,TO 2 ,A) 



(4tt) 2 



r/.r 



A(toi,to 2 , A) 



2 In 



A(toi, to 2 , to) 



(-A 4 - to 4 + 2m 2 h 2 )y 2 2y 3 (A 6 - 3to 2 A 4 + 3m 4 A 2 - to 6 ) ■ 



A 2 (TOi,TO2, A) 
d 4 q Z 4 



3A 3 (mi, TO2, A) 
/ A 2 - TO 2 \ 4 



(2^) 4 (p 2 - m 2 ){p 2 - TO 2 )(g 2 - to 2 ) \ q 2 - A 2 

1 



dx 



dy^6y{A 



to 2 ) In 



(4tt) 2 

-4y(A 2 -TO 2 ) + 6 A(mi,m2,A)lnf 
-A(toi, TO2, to) In '" 



A(mi, to 2 , to) 



A(toi,to 2 , A) 
1 

A(mi,m 2 , A)> 
(3y 2 A 4 - to 4 + 2to 2 A 2 
A(toi,to 2 , A) 



2/ 3 (A 6 -3to 2 A 4 + 3to 4 A 2 



A 2 (?7li,TO2,A) J 

q = l-p4X+p 3 (l -x-y), 

A(a, b, c) = a 2 (l — x — y) + 6 2 x + c 2 y + to 2 (.t 2 + y 2 — x — y + 2xy) + to 2 (.t 2 — x) 
+p z -p A {2x 2 + 2xy-2x). 



[1] BABAR Collaboration, B. Aubert et al, Phys. Rev. Lett. 87, 091801 (2001). 
[2] Belle Collaboration, K. Abe et al., Phys. Rev. Lett. 87, 091802 (2001). 

[3] The CDF Collaboration, F. Abe et al., Phys. Rev. D 58, 112004 (1998); Phys. Rev. Lett. 81, 2432 (1998). 
[4] M. Lusignoli and M. Masetti, Z. Phys. C 51, 549 (1991) 

[5] M.A. Ivanov, J.G. Korner and P. Santorelli, Phys. Rev. D 73, 054024 (2006). 
[6] M.A. Ivanov, J.G. Korner and O.N. Pakhomova, Phys. Lett. B 555, 189 (2003). 
[7] M.A. Ivanov, J.G. Korner and P. Santorelli, Phys. Rev. D 63, 074010 (2001). 

[8] P. Colangelo, G. Nardulli and N. Paver, Z. Phys. C 57, 43 (1993); E. Jenkins, M.E. Luke, A.V. Manohar, M.J. Savage, 
Nucl. Phys. B 390, 463 (1993); V.V. Kiselev, A.K. Likhoded and A.V. Tkabladze, Phys. Atom. Nucl. 56, 643 (1993); Yad. 
Fiz. 56, 128 (1993); V.V. Kiselev and A.V. Tkabladze, Phys. Rev. D 48, 5208 (1993); C.H. Chang and Y.Q. Chen, Phys. 
Rev. D 49, 3399 (1994); G.R. Lu, Y.D. Yang and H.B. Li, Phys. Lett. B 341, 391 (1995); Ikaros I.Y. Bigi, Phys. Lett. B 
371, 105 (1996); G.R. Lu, Y.D. Yang and H.B. Li, Phys.Rev. D 51, 2201 (1995); M.T. Choi and J.K. Kim, Phys. Rev. D 
53, 6670 (1996); C.H. Chang, J. P. Cheng and CD. Lu, Phys. Lett. B 425, 166 (1998); V.V. Kiselev, A.K. Likhoded and 
A.I. Onishchenko, Nucl. Phys. B 569, 473 (2000); AbdEl-Hady, A., J.H. Munoz and J.P Vary, Phys. Rev. D 62, 014019 
(2000); V.V. Kiselev, A.E. Kovalsky and A.K. Likhoded, Nucl. Phys. B 585, 353 (2000). 



10 



[9] V.V. Kisele v in N. Brambilla et al., CERN Yellow Report Heavy quarkonium physics, CERN-2005-005, Geneva: CERN, 
2005. 487 p. arXiv:hep-ph/0412158| C.H. Chang, Int. J. Mod. Phys. A 21, 777 (2006). 
[10] H.Y. Cheng, O.K. Chua and A. Soni, Phys. Rev. D71, 014030 (2005). 

[11] Z.Z. Xing and D.S. Du, Phys. Lett. B 270, 51 (1991); Phys. Lett. B 276, 511 (1992); D.S. Du, X.Q. Li, Z.T. Wei and B.S. 

Zou, Eur. Phys. J. A 4, 91 (1999); Y.S. Dai, D.S. Du, X.Q. Li, Z.T. Wei and B.S. Zou, Phys. Rev. D 60, 014014 (1999). 
[12] N.G. Deshpande, X.G. He, W.S. Hou and S. Pakvasa, Phys. Rev. Lett. 82, 2240 (1999); 

[13] N. Isgur, K. Maltman, J. Weinstein and T. Barnes, Phys. Rev. Lett. 64, 161 (1990); M.P. Locher, V.E. Markusin and 

H.Q. Zheng, Report No. PSI-PR-96-13 (unpublished); H. Lipkin, Nucl. Phys. B 244, 147 (1984); H. Lipkin, Phys. Lett. 

B 179, 278 (1986); H. Lipkin, Nucl. Phys. B 291, 720 (1987); H.J. Lipkin and B.S. Zou, Phys. Rev. D 53, 6693 (1996); 

P. Geiger and N. Isgur, Phys. Rev. Lett. 67, 1066 (1991); V.V. Anisovich, D.V. Bugg, A.V. Sarantsev and B.S. Zou, Phys. 

Rev. D 51, R4619 (1995); X.Q. Li, D.V. Bugg and B.S. Zou, Phys. Rev. D 55, 1421 (1997); X. Liu, B. Zhang and S.L. 

Zhu, Phys. Lett. B 645, (185) (2007). 
[14] M. Masetti, Phys. Lett. B 286, 160 (1992). 
[15] J.F. Liu and K.T. Chao, Phys. Rev. D 56, 4133 (1997). 
[16] Y.S. Dai and D.S. Du, Eur. Phys. J. C 9, 557 (1999). 
[17] V.V. Kiselev, J. Phys. G 30, 1445 (2004). 

[18] J.F. Cheng, D.S. Du and CD. Lii, Eur. Phys. J. C 45, 711 (2006). 

[19] M.A. Ivanov, J.G. Korner and P. Santorelli, Phys. Rev. D 73, 054024 (2006). 

[20] Y. Oh, T. Song and S.H. Lee, Phys. Rev. C 63, 034901 (2001) 

[21] M.A. Shifman, A.I. Vainshtein and V.I. Zakharov, Nucl. Phys. B 147, 385, 448 (1979). 
[22] X. Liu, X.Q. Zeng and X.Q. Li, Phys. Rev. D 74, 074003 (2006). 

[23] C. Itzykson and J.B. Zuber, Quantum Field Theory, McGraw-Hill Inc., New York, 1980. 

[24] M.E. Peskin and D.V. Schroeder, An Introduction to Quantum Field Theory, Addison- Wesley Publishing Company, New 
York, 1995. 

[25] W.M. Yao et al., Particle Data Group, J. Phys. G 33, 1 (2006). 

[26] The CLEO Collaboration, M. Artuso et al., Phys. Rev. Lett. 95, 251801 (2005). 

[27] C. Isola, M. Ladisa, G. Nardulli and P. Santorelli, Phys. Rev. D 68, 114001 (2003). 



